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Abstract 

The fits of 0ACa(NO 3 ) 2 0.6K(NO 3 ) (CKN) by schematic mode-coupling models [V. 
Krakoviack and C. Alba-Simionesco, J. Chem. Phys. 117, 2161-2171 (2002)] are analyzed 
by asymptotic expansions. The validity of both the power-law and the Cole-Cole-peak 
solutions for the critical spectrum are investigated. It is found that the critical spectrum 
derived from the fits is described by both expansions equally well when both expansions 
are carried out up to next-to-leading order. The expansions up to this order describe the 
data for 373K over two orders of magnitude in frequency. In this regime an effective power 
law of can be identified where the observed exponent a is smaller than its calculated value 
by about 15%; this finding can be explained by corrections to the leading-order terms in 
the asymptotic expansions. For higher temperatures, even smaller effective exponents are 
caused by a crossover to the alpha-peak spectrum. 
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1 Introduction 



For a given interaction potential mode- 
coupling theory (MCT) in its micro- 
scopic version can predict the transi- 
tion from a fluid to a glass [1]. Close 
to the transition, the glassy dynamics 
is ruled by universal scaling laws, that 
are independent of the details of the 
interactions. These universal laws have 
been applied frequently to describe 
experimental data, especially when 
the microscopic details are difficult 
to model, which is the case in most 
molecular glass formers. At larger dis- 
tances from the liquid-glass transition, 
non-universal corrections to scaling 



become important [2], These can be 
calculated exactly by asymptotic ex- 
pansion for a given MCT model, but 
introduce additional parameters for the 
fitting of data. An alternative to fitting 
higher-order terms of an asymptotic 
expansion is provided by schematic 
MCT models. These simplified models 
incorporate naturally both universal 
features and corrections, and can be 
used to fit data by a limited number of 
parameters. Similar to the microscopic 
models, asymptotic expansions can 
also be devised for the schematic mod- 
els for experimentally relevant param- 
eters, which in turn provides insight 
into the validity of the universal laws 



Preprint submitted to Journal of non-crystalline solids 



2 February 2008 



and the importance of the corrections 
for the data under investigation. 

One particular schematic model used 
frequently is the F\2 model which deals 
with a single correlator mimicking the 
density-correlators [3]. This model can 
be supplemented by a second corre- 
lator to describe some specific probe 
variable [4]; and this two-correlator 
model has been applied successfully 
to the description of experimental data 
[5,6,7,8,9,10]. Data for salol [11] and 
benzophenone (BZP) [12] were fit- 
ted recently by this model, and it was 
subsequently shown that an apparent 
violation of universal MCT formu- 
las [13,12] can be reconciled with 
a more detailed asymptotic analysis 
within MCT [14,15]: A Cole-Cole 
peak found earlier for this schematic 
model [16,17] was identified in the 
fits for these two substances. While 
the asymptotic approximation by the 
Cole-Cole peak improved the descrip- 
tion of salol only quantitatively, the 
presence of a Cole-Cole peak was es- 
sential for interpreting the BZP spectra 
even qualitatively. A remarkable wing 
in the Giga-Hertz regime could be ex- 
plained by a significant contribution of 
the Cole-Cole peak [15,18]. 

Another substance studied in detail is 
the molten salt 0. 4Ca (N0 3 ) 2 0.6K(NO 3 ) 
(CKN) which was measured in light- 
scattering experiments [19]. These data 
were successfully fitted by a number 
of different schematic models [20,21]. 
For these models the asymptotic ex- 
pansions will be derived in an effort to 
identify asymptotic features present in 
the data of CKN. The models will be 
introduced in Sec. 2, their asymptotic 
solutions presented in Sec. 3; Sec. 4 



discusses the two lowest temperatures 
considered in [21], Sec. 5 extends the 
analysis down to temperatures as low 
as in Ref. [20], and presents a conclu- 
sion. 



2 Schematic Model Fits 

MCT provides equations of motion for 
correlation functions § q (t) where the 
index q denotes the wave-number or a 
specific dynamical variable. With the 
initial conditions § q (t = 0) = 1 and 
d t ty q (t = 0) = 0, the MCT equations 
read 

+Q2 f dt'm q (t-t')d t ^ q (t > ) = 0, 

J 

(1) 

where Q. q and V q specify normal liq- 
uid dynamics. In the F\2 model, the 
dynamics of the density correlators is 
modeled by a single correlation func- 
tion with kernel 

m(t) =v l ty(t)+V2§(t) 2 . (2) 

The probing variable is given by a sec- 
ond correlator (|>a(0 and a correspond- 
ing kernel [4], 

m A (t) = v A ^(t)^ A (t). (3) 

The dynamics of ty q (t) is given by 
Eq. (1) with q replaced by A. A dif- 
ferent kernel for a second correlator 
was introduced by Alba-Simionesco 
and coworkers that shall be denoted 
by <|> r (0 [22], 

m r (t) = rm{t) = r [v^(t) +v 2 $(t) 2 ] . 

(4) 



2 



In addition, a number of effective cor- 
relators § s (t) were defined, some of 
whose were originally motivated to 
capture the DID mechanism [22,21], 

(M0 = 7*K0 2 +(1-Y)<M0 2 , (5a) 

k(0=YK0 2 +(i-Y)<kr(0, (5b) 

i(0 = ^(f) + (l-Y)M0 2 . (5c) 

i(0=#)+(i-T)M^ (5d) 

Here, <|>A,r(f) denotes a second corre- 
lator with the kernel from Eq. (3) or 
Eq. (4), respectively. The data was fit- 
ted with the susceptibility given by the 
Fourier transformation 

f'OO 

X" (to) = to / dt cos(cof)<MO • (6) 
Jo 

Further details and the values of 
the fit parameters can be found in 
Refs. [21,23]. The fit curves are re- 
produced in Fig. 1 together with the 
experimental data for CKN for the 
two lowest temperatures considered 
in [21]. Here, kernel (4) and the ef- 
fective correlator (5a) are applied; the 
frequencies are denoted in units of 
V = to/(27t). 



3 Asymptotic Solution for Schematic 
Models 



The general asymptotic formulas for 
the critical spectrum are given in the 
following and are specialized to the 
schematic models subsequently. From 
those formulas the asymptotic descrip- 
tion of CKN shall be inferred. The 
following discussion is restricted to 
spectra directly at the transition point, 
which will be called critical spectra. 



Such a critical spectrum is shared by 
larger and larger parts of the dynamics 
when the temperature is approaching 
the critical temperature T c . Two dif- 
ferent asymptotic approximations for 
the critical spectrum shall be used, a 
well-known power-law expansion and 
a modified Cole-Cole law introduced 
recently. 

The power-law solution for the critical 
spectrum is derived from the asymp- 
totic expansion for the correlators. It 
reads up to next-to-leading order [2], 

Xj(to) = [r(l -a)sm(na/2)]h q (m ) a 

x [l + k a K q {mo) a ] , 
(7a) 

with the prefactor 

k a = 2T{\-a) cos (na/2) /I, (7b) 

the critical amplitude h q , the critical ex- 
ponent a - given by the exponent pa- 
rameter X using the Euler-y function in 
I = T(l - a) 2 /T( 1 - 2a) - the micro- 
scopic time scale to, and a correction 
amplitude K q . The first line in Eq. (7a) 
constitutes the leading-order result: an 
increase of the critical spectrum with 
as (if. 

In contrast to Eq. (7), the modified 
Cole-Cole law is obtained from us- 
ing the expansion for the kernel m q (t) 
rather than the one for § q (t) [15]: 



+Z? (-/to/to<) 



2a 



(8a) 



where all parameters can be cast in a 
form using only results known from the 
power-law expansion. The amplitude is 
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given by 



(8b) 



with the glass -form factor f c q \ the char- 
acteristic frequency reads 



co^ =[(l-/, c )/(^r(l-a))] 



1/a 

(8c) 



and the correction amplitude is ex- 
pressed as 

K c q c =!+[(!- f q )/h q ]K q /X. (8d) 



For K" c 



0, the correction vanishes, 
and the leading-order result in Eq. (8a) 
is recovered as the original Cole-Cole 
function [27]. 

For the schematic models used in the 
following, one can rewrite all the pa- 
rameters needed in Eqs. (7) and (8) in 
terms of only two variables: the ex- 
ponent parameter X and the respective 
coupling parameters, v A or r. We ob- 
tain for the glass-form factors, 



fr = 



v A (l-X)-l 
va(I-X) 
r(l-X) 



r(l-X) + X 



(9a) 

for the critical amplitudes, 

h = X, h A =X/[v A (l-X) 2 }, 

rX (9b) 
[r(l-X)+X] 2 

and for the correction amplitudes, 
K=x, 



K A =K + 



l-X 

K r =K + X 2 



va(I-X) 
v A -\/(\-X) 
l-r 



-1 



r(l-X) + X' 



(9c) 



K can be determined straightforwardly 
from the exponent parameter X [17]. 
For the discussion below, K will al- 
ways contribute only a very small part 
to the correction amplitude; assuming 
K = would not alter any results sig- 
nificantly. 

The models for the second correlators 
in Eqs. (3) and (4) are in fact very sim- 
ilar, and the preceeding formulas also 
motivate a closer comparison. From 
setting f A = f r , one obtains a mapping 
between the coupling coefficients, 



v A X- 



i-r 



v A 



r 1 
X + Y^X' 



(10) 

When v A and r are related by Eq. (10), 
the ratio between the critical ampli- 
tudes becomes 



h r r(l-X) 



(ID 



which becomes unity for r = X 2 /(l — 
X) 2 . For that specific value of r, the 
mapping (10) is correct in leading 
asymptotic order: % A (u>) and 
have the same leading-order approxi- 
mation in both Eq. (7a) and (8a). Im- 
posing again Eq. (10), the ratio among 
the correction amplitudes is 



K a -k 
K r -K 



1 + 



r(l-X) 



(12) 



which cannot be unity with h a /h r si- 
multaneously. However, the correction 
amplitudes get closer to each other 
for larger values of r. Typical values 
for r and v A can be obtained from 
Ref. [21]; there, both models for the 
second correlator were used for fitting 
the CKN data independently. The val- 
ues reported in Ref. [21] agree within 
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a 15% margin with the mapping in 
Eq. (10); the values for r also allow 
for an approximate identification of 
the amplitudes, h A ~ h r , and the ra- 
tio in Eq. (12) lies between a maxi- 
mum of 2 and a minimum as low as 
1.2. Two conclusions can be drawn 
from the preceeding discussion. First, 
both models, %^((g)) and %"((&), are 
almost equivalent in their asymptotic 
properties for experimentally relevant 
parameter values; and second, the fit- 
ting procedures employed in Ref. [21] 
show this equivalence consistently. 

The asymptotic expansions for the ef- 
fective correlators and their spectra are 
calculated by inserting the results for 
§(t) and §A,r into Eq. (5) and retaining 
only terms up to second order. For the 
case (5a) one gets, 

fs = yf c2 + (i-y)fi r , 

h s = 2yfh + 2(\-y)f A , r h A ,r, 
K s =[y(h 2 + 2fhK) + (l-y) 

X {h\ r + 2f A ,rhA,r^A,r)\ /h s . 

The remaining cases from Eq. (5) can 
be calculated in the same way and yield 
similar formulas, so these results need 
not be given here. The parameters for 
the Cole-Cole law result from inserting 
the values from Eq. (13) into Eq. (8) 
with q = s. 



4 Results 

The results shown in this section will 
be restricted to temperatures at and 
above 383K as in [21]. There, this re- 
striction was motivated to avoid the 
artefacts in the original data [19] identi- 
fied for lower temperatures [24,25,26]. 
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Fig. 1. Power-law expansion for CKN. 
The full curves reproduce the data from 
Ref. [19]. The dashed-dotted curves 
present schematic-model fits using ker- 
nel (4) and the effective correlator (5 a) 
with r = 8.6 as in [21]. The full lines 
labeled c display the critical spectra for 
X = 0.656 (a = 0.344, upper panel) and 
X = 0.7 (a = 0.327, lower panel), respec- 
tively. The dotted lines labeled of and 
the dashed lines labeled (s) 2a show the 
leading and next-to-leading order result 
of Eq. (7a). The points where the critical 
spectra deviate by 10% from the approxi- 
mation by (7a) are indicated by circles. 

However, the results derived for the 
critical spectrum in the present section 
will be shown to be relevant also for 
lower temperatures in Sec. 5. For the 
discussion of the critical spectrum, an 
appropriate point on the transition sur- 
face of the F\i model needs to fixed. In 
the F\2 model such a critical point can 
be specified uniquely by the value of X. 
One reasonable choice is a transition 
point close to the fit parameters for the 
lowest temperature under discussion. 
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Such a point is chosen for the critical 
spectrum shown in the upper panel of 
Fig. 1; the exponent parameter for this 
critical point is X = 0.656. Another 
way to determine a suitable transition 
point is suggested by an extrapolation 
procedure in Ref. [21], where X ~ 0.7 
is proposed for the description of the 
dynamics. The critical spectrum for 
the transition point satisfying X = 0.7 
is added to the lower panel of Fig. 1; 
to show also experimental results close 
to that transition point, data and fit of 
the nearby 393K spectrum is added. 

The asymptotic approximation by 
Eq. (7a) is shown in Fig. 1, and it 
is seen clearly that the leading-order 
power law describes the critical spec- 
trum only for v < 1GHz, regardless of 
the chosen transition point. Including 
the correction describes the critical so- 
lutions up to 300GHz and 150GHz, re- 
spectively. However, Eq. (7a) explains 
the measured spectrum for 383K for 
only about one order of magnitude 
in frequency. The apparent power-law 
behavior v a with a' = 0.22 between 
2GHz and 200GHz for 383K does in 
turn not show the asymptotic value 
a = 0.344 for the critical exponent. 
This is because the separation of this 
state point from the transition is too 
large. 

The alternative asymptotic approxima- 
tion by the Cole-Cole law, Eq. (8a), 
is shown in Figure 2. Similar to the 
power-law solution, the leading-order 
Cole-Cole formula applies to the crit- 
ical spectra only for v < 1GHz, but 
has no relevance for the description of 
the present data. The formula includ- 
ing the correction describes the critical 
spectra in upper and lower panel for 
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Fig. 2. Cole-Cole law for CKN. The full 
curves display the experimental data [19] 
as before. Also reproduced from Fig. 1 are 
the respective critical spectra (c) and the 
power-law approximations by Eq. (7) up 
to next-to-leading order (co 2<3 ). The dotted 
line labeled cc and the dashed line labeled 
ccc show the leading- and next-to-leading 
order approximation by the Cole-Cole law, 
Eq. (8a) for frequencies indicated by the 
arrows. A 10% deviation of the critical 
spectra from their respective approxima- 
tion by Eq. (8a) are shown by the triangles. 

v < 200GHz and v < 150GHz, respec- 
tively. The comparison of both panels 
also reveals that the choice of a par- 
ticular transition point does not lead 
to major differences. In both cases, 
Cole-Cole and power-law expansion 
are practicably indistinguishable for 
v < 500GHz. This is also the regime 
where they both approximate the criti- 
cal spectrum well. 

As seen in Fig. 1, the range of va- 
lidity for approximation (7a) increases 
by two decades when the corrections 
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2750 
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v cc 
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0.496 


0.491 


0.483 


0.478 


0.510 


0.507 


k c s c 


0.63 


0.68 


0.67 


0.68 


0.67 


0.59 



Table 1 



Selected parameter values for the asymptotic approximation of the critical spectra. The 
designation of the cases I-VI is the same as in [21]; the first line shows the definition of 
the effective correlator § s . 

are included, K s = —0.31 for the up- 
per panel; such a value for K s can be 
considered moderate compared to BZP 
and salol where the corresponding val- 
ues are around five times larger [15]. 
For the lower panel K s = —0.44. On 
the contrary in Fig. 2, the corrections in 
(8d) for approximation (8a) are much 
larger for CKN, £f = 0.63 and K c s c = 
0.54 for upper and lower panel, than 
for BZP and salol where the corrections 
are always close to zero for the Cole- 
Cole law. Table 1 lists the values of 
characteristic parameters for the mod- 
els used in Ref. [21]. Model I is the one 
discussed in detail above for Figs. 1 
and 2. Most parameters are remark- 
ably similar throughout all the mod- 
els. The location of the Cole-Cole peak 
V c s varies by only 30%, so the scenario 
shown in Fig. 2 is typical for all these 
models. The values for the case X = 
0.7, corresponding to the lower panels 
in Figs. 1 and 2, are very similar and 
are therefore omitted. 
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Fig. 3. Spectra for 373K and 363K. Crit- 
ical spectrum (c), power-law approxima- 
tion (co 2a ) and Cole-Cole law (ccc) are re- 
produced from Figs. 1 and 2. In the up- 
per panel, the squares indicate a 10% de- 
viation of the approximation (7) from the 
data. An effective power law CO 3 is shown 
as dotted line. 

5 Discussion 



Figures 1 and 2 indicate that both the 
approximations, Eq. (7a) and (8a), as 
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well as the full critical spectra do not 
describe significant parts of the mea- 
sured data for T ^ 383K. For lower 
temperatures, two additional spectra 
with fits from Ref. [20] are available, 
T = 373K and T = 363K. The data for 
these temperatures from Ref. [19] are 
shown in Fig. 3. It is apparent from 
the upper panel of Fig. 3 that the spec- 
trum for T = 373K coincides with the 
critical spectrum for 2.5 decades; the 
asymptotic expansions describe the 
data for over two orders of magnitude. 

The corrections to the power law, cf. 
Eq. (7 a), lead to a renormalization of 
the asymptotic power law ro 344 to an 
effective power law (O 3 . In contrast 
to the findings for salol and BZR the 
Cole-Cole formula (8a) is not essential 
for understanding these smaller expo- 
nents for CKN [14,15]. The asymptotic 
expansions also reveal that the mini- 
mum seen in the data for CKN can 
be regarded as a [3-minimum with a 
von Schweidler wing for lower and a 
critical spectrum for higher frequen- 
cies. Therefore, for temperatures T > 
383K, earlier analyses of the minimum 
can be justified [28]. 

For the 373K spectrum a fit can be 
found with the same parameters as for 
T = 383K by adjusting vi and V2 only. 
The resulting state is in fact close to 
the critical point for X = 0.656, and a 
similar point was found earlier for the 
same temperature for slightly differ- 
ent parameters [20]. The situation for 
T = 363K is different. The spectrum 
shown in the lower panel of Fig. 3 devi- 
ates drastically from the critical spec- 
trum for v < ITHz. Hence, in this case 
also the approximations to the critical 
spectrum are not applicable. Although 



the data for T = 363K could be fitted 
in [20], this was possible only by mov- 
ing away from the path for higher tem- 
peratures rather considerably [21]. One 
might interpret that as an indication of 
the mentioned experimental artefacts 
[21,24,25,26]. 

For temperatures T ^ 363K, the expo- 
nents reported in Fig. 4 of Ref. [25] 
are all equal or larger than 0.4. These 
temperatures, when interpreted within 
MCT, are therefore located below T c . 
This is a consequence of the evolu- 
tion of the spectra as seen for higher 
temperatures, as well as the asymp- 
totic analysis at the critical point: Both 
mechanisms only allow for smaller 
effective exponents compared to the 
asymptotic value for a which cannot 
be larger than 0.396. Below T c , larger 
effective exponents are not in contra- 
diction to MCT; however, they do not 
describe a critical spectrum. 

In conclusion, for T ^ 363K, in the 
spectra from [19] as well as [25] no 
critical MCT spectrum can be demon- 
strated; T = 383K marks the highest 
temperature where part of the critical 
spectrum is observable; and the data 
for 373K, cf. upper panel of Fig. 3, ex- 
hibit the critical MCT spectrum rather 
convincingly. 
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